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ABSTRACT 


Experimental tables of availabilities at time t are 
obtained for a device whose performance is described by an 
alternating renewal process with a finite number of failure- 
renewal cycles, corresponding to having a finite spares 
backlog. Failure and repair rates are assumed to be constant, 
and attention is restricted to cases in which the repair rate 


1s larger than the failure rate. 
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I. INTRODUCTION 


The most commonly encountered working definition of the 


Mavailability”™ of a device, 


MTTF 
Gl) Availability = ————— 
MTITF + MTTR 
represents the long-term or steady-state probability that 
the device will be found in an "up" or functioning condition 
when two specific conditions are satisfied. One condition is 
that there is an alternation of failure and repair cycles 
in which times to failure and times to repair are independent 
realizations from some failure and repair distributions 
satisfying minimal regularity conditions. The second, and 
here Most important condition, is that the alternation of 
failure and repair continues indefinitely, so that the 
performance of the device is described by a standard alternating 
renewal process. 

For Many equipments, the second condition cited above 
implies access to an infinite backlog of spares. In many 
operational contexts this sort of spares support cannot be 
realistically assumed. 

If availability is considered to derive from an alter- 
nating renewal process with a finite number of cycles, 
corresponding to a finite backlog of spares, then expressions 


for availability become complex as compared with equation (1.1). 





This thesis is devoted to computational experiments with 
some "finite eaeracn availability expressions. The end 
objective of such experiments is to be able to determine the 
circumstances in which equation (1.1) furnishes an adequate 
approximation, or alternately to be able to provide compu- 


tationally feasible alternatives to its use. 





Il. MATHEMATICAL MODEL 


The mathematical model on which the usual expressions 
for the availability of a device are based is an alternating 
renewal process with an infinite number of failure-repair 
Cycles. In situations where repair requires replacement of 
a failured device by a spare, this corresponds to having an 
infinite number of spares. The model studied here is modi- 
fied to allow only a finite number of failure-repair cycles, 
corresponding to having a finite number of spares. 

The simplest assumptions about failure and repair times 
are made; failure rates are constant, and repair rates are 
constant. Only those processes that begin with a functioning 
device installed are considered. 

In greater detail, the failure-repair process considered 


1s as shown in figure 2.1, 





IS S. | end 


ee RS a 


Peueres cael. Failure-repair process. 





where n is the number of spares, Ty is the time to failure 
for the original device, Tl; T i are the times to 
failure for the n spares, and Sq: Sy» re Snel are the 
times to replace the original device and the first n-l spares. 
It is assumed that T), Sys ceeeee Tel?  n-1? TL are inde- 
pendent random variables, and that Tosseeees , T, are 
exponentially distributed with failure rate ’, while 

SQ» \ oer Sn-] are exponentially distributed with repair 
rate n. 

The availability at time t, A(t) of the original 
device, supported by its backlog of n spares, is the | 
probability that the process shown in Figure 2.1 is in an 
"up" condition at time t; i.e., that at time t either the 
Original device or one of its spares is installed and still 
functioning. 

The increment in availability at time t due to the 
eee Spare, Ty (t) is defined by 
2.1) Tee) = ee) ae Ay _ 4 (t) | Sane arene PE Oe 

Before proceeding to a derivation of expressions for 
I(t) and A(t) in a general case, two special cases are 
considered; repair rate n equal to infinity, and repair 
rate n equal to failure rate »}. These are boundary cases 
for the cases of likely practical interest, in which it is 


reasonable to expect that repair rate will exceed failure 


na be , 





In any case, Ay(t), availability at time t with no 


Spares is given by 


(2.2) Aj(t) = PIT, > th = e**, t > 0. 


In the following sections, it will be convenient to let 


(2.3) Gt ye Badaee + T, 
2 Sy ean ons 
eo eG 
= (S, *T)) + ....-- (See tem Ta) 


A. REPAIR RATE EQUAL TO INFINITY 
The simplest case is the one in which no time is required 
to repair a failed unit, provided a spare unit is available. 


In this case the contribution of the first spare is 
c 
(2.4) T(t) nf P(T, > t-s|U, = $s] ch 
0 


where 


io s > 0 


3 


f.. (s) = Ae 
Uy 


1s the gamma {1,\} density. Thus 


ie 
12.5) T(t) | ee ea ds 
0 





_ re AS t, 


and the availability of a system having one spare is 
(2. 5) A, (t) = Ay (t) + I, (t) 
= 5 al + Xt) 


The contribution of the second spare is 


iG 
i) I, (t) is P[T, > t-s|U, = s)%y, (8) ds , 
0 


where 


ae a) 8 


is the gamma {2,A} density. Thus 


te 
Z__-AS 
_ -A(t-s) A“se 
(2.8) (ae) [ e —ty- 4s 
0 
Z.-At L 
= re) Cc VAG : 


and the availability of a system having two spares is 
2 
(2.9) ies “Oe ape UE 


: h 
Generally, the contribution of the n* spare 1s 


10 





s 
© 300) I(t) = Bee > t-s|U__4 = sjfy (yds 
Q. : 
where 
yh n-li-As 


£ (5.) S715 60 eee ei 
ened 7 


is the gamma {n,A} density. Thus 


c 
(@.11) it G3) oe eta Ae 
0 Penn 


Heee = A el 


and the availability of a system having n spares is 


2), A(t) = eo “PCA + at + 


oe @ @ # @ 


Cole 
ih as 


Availabilities A(t) Obtained from equation (2.12) are 


Shown in table 1. In this table n represents the number of 


Spares, and the contribution I(t) of the ace 


found by subtracting A,-1 S*) from A(t). The availabilities 


Shown in the last column are for an alternating renewal pro- 


(At 
Ti 


3 


a 


spare can be 


cess with an infinite number of failure-repair cycles. 


corresponds to having an infinite number of spares. 


B. REPAIR RATE EQUAL TO FAILURE RATE 


The failure-repair process considered is that in which 


failure times and repair times are exponentially distributed 


With equal rates, i.e., A =n. 


ial 
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In this case, the contribution of the first spare is 


=c 
(2.13) 5 (C2) -| P[T, > t-s|W, = s]£y_ (8) ds , 
0 


where 


ee er PS 


is the gamma {2,\} density. Thus 


ie 
a Geass Neca ds 
0 


12—_7At Z 


t 


(14) I, (t) 


and the availability of a system having one spare is 


Be)! AGO) = LaNCt) 


2 
ot + Gey") 


th 


Generally, the contribution of the n Spare is 


16 
fe... 15) I(t) f es > t-s[W_ = s]fy (s) ds , 
0 n-1 


where 
,2n,én-1_-As 


18 (en) = Smee a ose 2 es ee 
ea | n ~ 
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is the gamma {2n,A} density. Thus 


i -AS 
(2216) ae) f Py GE mS)) Tomy (as) 7 re: 
0 


1B 
2n -At 
Ne [ -fn-1 Ae 
2n- 
0 


yon,“ 


(2n)! ? 


and the availability of a system having n spares is 


_ 2 2n 
(2.17) A(t) =e MOT Ca ie : i (At) 


C. REPAIR RATE GREATER THAN FAILURE RATE 

The failure-repair process considered in this section 
1S visualized as one in which the repair rate is greater 
than the failure rate. This influences the format in which 


the results are displayed. 


In this case, the contribution of the nth Spare 1s 


© 
(22 3) T(t) “f P[T. > t-s|Wo oy = s] fy (s) ds , 
0 n-l 


where 


ie 
£ (so) = ti Cass (t-s) ds , 
ae [ Coed ea 


14 





000S° 000s" 
000S° 000s" 
000S° 000s" 
z00S° zoos" 
ztos*  ztos" 
Z60S° 2605" 
4005 YESS 
9TI9°  9TT9" 
6589° 689° 


cl08- UNS: 


ayy G@qllky 


000S° 
000S° 
000S © 
CO0O0S © 
CTOS° 
C60S° 
ELIS: 
9TT9° 
6£89° 


cl08° 


(2) °y 


L66t- 
0000S" 
0000S" 
COOS © 
Guns: 
C60S © 
ELS: 
9TT9° 
ete 


cl08- 


(2) ®y 


£86t 
9667 © 
000S° 
CUUS. 
cU0S © 
C60S ° 
CLO 
9TT9" 
6£89° 


cf08- 


eney 


CL6V~ 
VL6V" 
S660 - 
TOOS~ 
CCOS © 
€60S " 
LLOS° 
9TT9° 
6£89° 


co08° 


(2) ?v 


8r9Or- 
T980- 
T967- 
S660 - 
CLUS 
C60S° 
LEO Se 
9LT9° 
6£89° 


cl08- 


(1) °v 


6S65° S¢E9z° SbrITt’ ¢&Z0° L 
Stryp° OTS” 608T° TLPrO’ 9 
6LZv° 0 Oo L@Tb.OCOU 997" COO T60° S 
Cv6v° vv9or” 209c° B8r9oTt’ v 
vOOS” ¢dé6bv° 6Tttv’ S8ElZ2° c 
T60S° ¢80S° 2d96r° O90F* C 
ZZ9S° ZLL9S° TL9S° 8gTss° I 
9TT9° OTT9° PTT9"* 2S09° SZ 
689° 6£89° 6£€89° ¢289° 0S° 
cfos° €¢F80° ¢ClO8° TE08s° SZ- 
CO Oe ey ents ry 
(2) “v 
v c C T U 


u 
1¥ = 2U osed oy 20OF (4) V Z OTGRL 


15 





and 


is the gamma {n,n} density, 


while ’ 


1g 
Ey ne 7 — a eee 


h- 


1s the gamma {n,A} density. Thus 


ic S 
Tl | 7 7 
0 0 


: Ces) du ds - 


Inverting the order of integration, equation (2.19) 


becomes 
iG t 
n n 
-A (t- -j] - r -] 
(2 20) ft) -[ f e (t-s) Tay we e es Tum) fous 
QOYu 


. 5 ND. ds du 


IE : ie 
ile, 18 
= eee pit. (M-Aju isa = ds du. 
ON Gen a: 0 ‘i 


Let v= s - u. Then 


1E t-u = 
| (s-u)2 + ds a) yn-l ay = ae 
0 0 


16 





Thus equation (2.20) reduces to 


1e 
n,n 
(2.21) I(t) =e “a fu Ete Pe 14 au 
0 


and the availability of a system having n Spares is 


222) A, (t) = Ay (t) ft [,(t) a oe ab ri (t) 


Note than when n = A, equation (2.21) reduces to 


equation (2./6), since then 


At ne n 
0 


Let u= tv. Then 
t 
wrt ee uy? au 


ae 
{ ae a tdv 
0 


t 
sal vet qev)™ av 
0 


z ,2n Gr) sent 1) 
NC2ne 
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so that 


a tay 
I : A, 2n Pn) TP (nt) 
n {t) n nt nt 
ie) ea 
n): 
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III. APPROXIMATION OF THE MATHEMATICAL MODEL 


In section II, we derived a mathematical model for 
aVailability. 
In this section, we discuss methods of approximation to 


obtain numerical values of availability. 


A. EXPONENTIAL EXPANSION APPROXIMATION 


The integral in equation (2.21) can be approximated by 


expanding its exponential term, i.e 


eae 


2 
ea) ene en-us + Ss a 


Thus the integral becomes 


t 
Z 
(3.2) f uP *ct-u)™ {a-(n-ayet uw? . } du 
0 
c 
a Tee eune du 
0 
t 
eee) ecu 
0 


We know that, 


19 





2 


ut-t eu Ae ,2n Beye! 
n 
0 
Thus equation (3.2) becomes 
-At_n,n 
_ @ na an Pin) Piney) 
(3.5) Hi Sted T (n)P(Mmel) {t '(2n+1) 


OTe erica ea tice) 
ee ene = 


. (n-n)207B?4 Pentz)P (ntl) 


spine, 


ee n 19 f 
ont) At) ia) nea) Eeath) 


. (nt-At)? P(nt2) 
ae. 


nr 


Computational experiments, with the approximation 
represented by equation (3.3) have indicated unsatisfactory 
convergence behavior when (nt-At) is large, a case of some 


practical interest, and so this approach was not pursued. 


B. SIMPSON'S RULE APPROXIMATION 
The integral in equation (2.21) can be approximated by 


using Simpson's rule. 


20 





(3.4) Hoe (eeume 0 ae du 


n-1 
(tV) {t(1-v)}" pracile ee dv 
0 


uk 
,2n po Sian e7 (Nt-At)v eee. 


Thus equation (2.21) becomes 
n n t 
nt Kt n n -(nt-at)v 
(a5 ) I(t) = @ At eptaes f vole) -¢ dv 
0 


Now Simpson's rule can be applied to the integral in 
equation (3.5) to obtain numerical values of availability. 


Simpson's rule as applied is 


b 
aig 
5x20 } £(x)dx = RiGee gt a ey ea Tey ty a) 
a 
where , « ba 
m 3 
and Xi = 4, X, Seinecact Ilene pate ot othe ; Xatl = atmh = b, 


Blt 








while 


aia 


£(X,)5 ¥2 * 
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IV. TABLE 





In this section experimental tables of availabilities 
A(t) are shown, which were obtained from the mathematical 
model evaluated using Simpson's rule, i.e., using equation (3.6) 
to evaluate equation (3.5) with m = 500, X, = 0.0001 and 
Xe] =0.9999. The reason for choosing m = 500 is that compu- 
tational experiments with choices of m greater than or equal 
to 500 gave a stable and accurate result. 

In these tables n represents the number of spares, 


and the contribution of the nth 


spare I(t) can be found by 
subtracting A. 6t) from A(t). 

The availabilities shown in this section, Table 3-12, 
are for cases in Which nt > At. 

The availabilities shown in the last column are for an 
alternating renewal process with an infinite number of failure- 


repair cycles. This corresponds to having an infinite number 


of spares. 
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V. SUMMARY AND CONCLUSIONS 


Certain computational approaches were tried for obtain- 
ing availabilities for a device supported by only a finite 
backlog of spares, using the simple assumptions that failure 
and repair rates are constant. Real failure and repair distri- 
butions may be more complex, but the case considered is a 
good case for initial computational experiments. 

Qf the two approaches tried, neither proved entirely 
Satisfactory in obtaining availabilities in a way that is 
fast and suitable for use with small-scale computational 
facilities, e.g.. hand-held calculators. Also neither was 
effective. over ane entire range of failure rate and repair 
rate combinations that might be of interest. 

Since an easily used, readily accessible,way to assess 
the impact of finite spares backlogs on availability is 
desirable in many mission planning contexts, further computa- 
tional approaches should be tried. 

The tables presented in section IV give availability 
values with which the result of such experiments can be 


compared. 
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